Thermodynamic modeling of the LiF-YF3 phase diagram 



I. A. dos Santos 3 *, D. Klimm b , S. L. Baldochi 3 , 1. M. Ranieri 3 

"Instituto de Pesquisas Energeticas e Nucleares, CP 11049, Butanta 05422-970, Sao Paulo, SP, Brazil 
b Leibniz Institute for Crystal Growth, Max-Born-Strafie 2, 12489 Berlin, Germany 



Abstract 

A thermodynamic optimization of the L1F-YF3 binary phase diagram was performed by fitting the Gibbs energy 
functions to experimental data that were taken from the literature, as well as from own thermoanalytic measurements 
(DTA and DSC) on HF-treated samples. The Gibbs energy functions for the end member compounds were taken from 
the literature. Excess energy terms, which describe the effect of interaction between the two fluoride compounds in 
the liquid phase, were expressed by the Redlich-Kister polynomial function. The calculated phase diagram and ther- 
modynamic properties for the unique formed compound, LiYF4, are in reasonable agreement with the experimental 
data. 
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1. Introduction 

The main interest of the system LiF- YF3 is due to the 
technological applications of L1YF4 crystals (YLF) as 
active laser medium Thomaef al. determined the 

phase diagram of this system using differential thermal 
analysis (DTA) data. A eutectic reaction was observed 
for the composition of 80mol% LiF-20mol% YF3 at 
963 K and a peritectic reaction for the composition of 
49mol% LiF-51 mol% YF 3 at 1098K fl. 

Harris et al. |4j performed an extensive study of the 
melting behavior of YLF by DTA and microstructural 
studies, it was assumed that YLF has a congruent melt- 
ing behavior. Nevertheless in the growth of crystals by 
the Czochralski method, YLF was easily seeded only 
with the composition of 50.7 mol%LiF-49.3 mol% YF3. 
They concluded that the degree of congruence of this 
compound is closely related to contamination with oxy- 
gen compounds, mainly with moisture J^H^]. Then it 
was proposed a phase diagram of the system LiF-YF3 
composed by two eutectics, one with composition of 
80mol% LiF-20mol% YF 3 at 979(2) K, and the other 
to the composition 49 mol% LiF-5 1 mol% YF3 at tem- 
perature 1 103(2) K, the same temperature was attributed 
to the melting of YLF. 

Although the experimental phase diagram of this sys- 
tem is well known, there is scarce data concerning ther- 
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modynamic properties of fluorides. Regarding the rare 
earth trifluorides, including YF3, Barin [7] collected the 
commonly used data. Lyapunov et al. obtained YLF en- 
thalpies and heat capacities in the solid and liquid state 
by the mixing method using a massive calorimeter |0]. 

In this work, the L1F-YF3 phase diagram was opti- 
mized by the fitting of the Gibbs energy function, taking 
into account both the experimental data from the liter- 
ature 10] and the data obtained in this work using the 
differential thermal analysis (DTA). Gibbs excess en- 
ergy terms in the liquid phase, which describe the ef- 
fects of interaction between the two fluorides, were ex- 
pressed by the Redlich-Kister polynomial function 
Heat capacity (C p ), enthalpy of formation and entropy 
at 298.15 K were assessed for the YLF compound and 
compared with those reported in the literature. 



2. Experimental 

Samples were prepared using commercial LiF (Aid- 
rich, 99.9%) purified by the zone melting method, and 
YF3 synthesized from the oxide (Y2O3, Alfa Aesar, 
99.99%) by the hydrofluorination method [10], under 
reactive atmosphere of HF and Ar in both cases. DTA 
curves were obtained using a TGA-DTA equipment 
from TA Instruments, model 2960. The experiments 
were performed under Ar flow, using Pt/Au crucibles; 
heating rate of lOK/min and samples masses around 
50 mg. The melting point of the pure substances, the 
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temperature of the solid phase transformation and of the 
invariant reactions were considered as the extrapolated 
onset of the thermal event. The liquidus temperatures of 
the intermediary compositions were evaluated from the 
extrapolated offset temperatures. 

The heat capacity of the YLF was evaluated using a 
Netzsch STA 409 PC Luxx heat-flux differential scan- 
ning calorimeter (DSC), where a proper Cp carrier was 
installed. DSC setup was temperature calibrated mea- 
suring the melting points of Zn, Au, In, Ni, and the 
phase transformation of BaCO^. All experiments were 
carried out under Ar flow of 50 cm 3 /min and the sam- 
ples were place in Pt/Au crucibles with lid. YLF heat 
capacity was determined according to ASTM-E-1269 
method. In this method, the sample DSC heat flow sig- 
nal is compared to the DSC signal of a calibration stan- 
dard of known specific heat (sapphire in this case). Both 
curves are corrected by a baseline correction experiment 
where empty reference and sample crucibles are placed 
in the DSC furnace and the system signal drift is mea- 
sured under identical experimental conditions. The ex- 
periments consisted of three steps, an isothermal seg- 
ment in 40 °C for 20 min, a dynamic heating segment 
with a lOK/min heating rate and a final isothermal seg- 
ment at maximum temperature for 5 min. 

3. Thermodynamic method 

To describe a T - X binary phase diagram it is nec- 
essary to define the Gibbs energy functions for all com- 
pounds in the system and the Gibbs functions of mixing, 
if solution phases are present. Usually these functions 
are unknown for most of the solutions, thus a thermody- 
namic assessment is required in order to determine the 
excess energy of the mixing. The Gibbs equation for the 
compounds is defined as function of enthalpy and en- 
tropy at the reference temperature state (298. 15 K) and 
can be obtained from the heat capacity function (Cp(T)) 
as follows: 



G(T) = H(T)-S(T)T 

H(T) = H° mA5K + f C P dT 

J298.15K 1 



(1) 

(2) 
(3) 



The Cp(T) functions can be obtained by fitting a set 
of experimental data at a suitable polynomial function 
expressed by: 



Depending on the compound, more terms in this 
Cp(T) expression are added or disregarded in order to 
obtain the best fitting. The thermodynamic data for the 
end members LiF and YF3 were taken from the com- 
pilation by Barin 10]. According to equation ([TJ and 
considering the equations © and (0, the enthalpy of 
formation, the absolute entropy at the reference temper- 
ature and the heat capacity are required to determine 
the minimum of G and thus thermodynamic equilib- 
rium. For the LiYF4 intermediate compound these data 
are not available, therefore the Neumann-Kopp rule was 
assumed to set the initial Cp equation and calorimetric 
data, afterwards these latter parameters were properly 
assessed by optimization. 

For the liquid solution phase, function G(T) is ex- 
pressed as the sum of the Gibbs energy weighed con- 
tribution of the pure compounds (Go), the contribution 
of an ideal mixture (G/o) and finally a term related to 
the non-ideal interaction, defined as the excess energy 
(G ex ). The sub-regular solution model of Redlich-Kister 
was adopted to describe the excess energy of the liquid 
phase in this system |U. For binary systems, the poly- 
nomial model is expressed by: 



G e x = x A x B 2^ Lj (x A - x B ) J 



(5) 



Cp(T) = a + bT + cT 



-2 



(4) 



where x A and xg are the molar fractions of components 
A and B , respectively. Lj terms represent the interac- 
tion coefficients between the basis compounds and they 
are given as a linear function of temperature. The opti- 
mization was performed using the OptiSage module in 
the FactSage 6.2 software [12], which uses the Bayesian 
Algorithm [13]. This algorithm is based on a probabil- 
ity model to obtain the fit between the theoretical Gibbs 
energy functions and the experimental data. 

4. Results and discussion 

4.1. Experimental results 

DTA curves for pure LiF, for three mixed LiF- 
YF3 samples with different composition, and for pure 
YF3 are shown in Figure Q] The curves for LiF and 
LiF+20 mol% YF3 exhibit only one endothermic peak 
with onsets near 1115 or 975 K, respectively. These 
peaks are due to LiF melting (« 1115K) and the 
LiF/LiYF4 eutectic 975 K, the eutectic point is close 
to this composition, cf. Figure [3). DTA curves for 
35 and 60 mol% YF3 show two endothermic events, the 
first peak characterizes an invariant reaction (e.g. eutec- 
tic or peritectic reaction) and the second broader peak 
marks the end of melting of the primary phase (YLF and 
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Table 1: Aff(298.15K) (kJmor 1 ), S(298.15K) (JKT 1 mol 


), AHf (kjmol ) and Cp data for LiF, YF3 and the intermediate compound LiYF4. 


Compound 


AH(298.15K) 


S (298.15 K) 


AH /(1) I 


AH f(2)L 


a 


b 


c 


LiF(S) 


-616.931 


35.660 


27.09 [7] 


27.68 


42.689 


1.742 x 10" 2 


-5.301 x 10 s 


LiF(Z) 


-594.581 


42.997 






64.183 






YF3GSI) 


-1718.368 


109.960 






99.411 


7.427 x 10~ 3 


-5.690 x 10 5 


YF 3 (S2) 


-2161.269 


-960.528^ 


27.97 [7] 


29.79 [11] 


-319.448 


2.128 x 10- 1 


2.818 x 10 8 


LiYF^ 


-2355.780 


138.325 


67.65 [8] 


63.538 


142.101 


2.484 x 10" 2 


-1.099 x 10 6 



"Data taken from Barin 01 and from Lyapunov et al. (g). 
'DSC data measured in this work and from fTlll . 
^Extrapolated value for the high temperature phase. 

f, A//(298.15 K) and 5(298.15 K) were assessed in this work and Cp function was estimated by the Neumann-Kopp rule. 



YF3 respectively) at the liquidus. In the YF3 curve the 
solid state transformation (1338 K) and fusion (1403 K) 
can be recognized very clearly. These two thermal 
events show very similar peaks in area, therefore sim- 
ilar enthalpy is required in these transformations, as it 
was observed previously [11]. 
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Figure 1: DTA curves for several LiF-YF^ compositions. 

The heat of fusion AHf for LiF, YF3 and OYF4 was 
calculated from the melting peak area of these com- 
pounds. Reasonable agreement with the literature was 
found in the AHf (Table[T]i. Nevertheless the Neumann- 
Kopp rule had been considered to set YLF Cp(T) func- 
tion used to perform the assessment, Cp data for this in- 
termediate compound were also measured by DSC. Fig- 
ure [2]compares DSC measured data and those estimated 
by the Neumann-Kopp rule. One can see that for the 
considered temperature range, there are no strong devi- 
ations (not bigger than 3%) between the experimental 
and estimated Cp data. Moreover, taking into account 
the polynomial fitting of the DSC experimental data us- 
ing the function (3J, the a and b parameters obtained are 



very similar to those achieved by Lyapunov et al. 18] 
and the difference is less than 3% for the c parameter. 
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Model: C p (T)-a+bT+c/T 
Chi A 2/DoF = 0.26963 
R A 2 = 0.99196 
a= 161±2 
b = (9±2)E-3 
c = (-3.6±0.2)E+6 




□ C p data obtained experimentally by DSC 
o C p data estimated by Neumann-Kopp rule 
Fitting considering C„-a+bT+c/T 2 model 
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Figure 2: Cp data measured experimentally through DSC technique 
and Cp estimated by Neumann-Kopp rule. 



4.2. Thermodynamic assessment 

A proper optimization of binary systems using a poly- 
nomial model to represent the Gibbs energy for the solu- 
tion phases is strongly dependent on quantity and accu- 
rate of experimental data available for the system stud- 
ied. Therefore, the experimental phase diagram of the 
system LiF-YF3 was constructed taking into account 
the experimental data evaluated by Thoma et al. ^ 
(squares), and the DTA data obtained in this work (stars) 
(Figure 0). It can be seen that our experimental points 
are in agreement with those collected from literature. 

The LiF-YF3 phase diagram has been optimized ac- 
cording to the Redlich-Kister polynomial model using 
the Bayesian Optimization Algorithm of FactSage lfl2ll . 
The data available in the literature for this system are 
data of solidus and liquidus lines (experimental T - X 
phase diagram) and Cp and calorimetric properties for 
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the end members (LiF and YF3). No enthalpies of mix- 
ing or activity data were considered for the LiF-YF3 
liquid phase. Consequently many degrees of freedom 
are present in the assessment, making the setting of the 
appropriately solution features not straightforward. It 
was observed by other authors working with different 
LiF-L«F3 systems (Ln = La-Sm; these systems do not 
contain a LiLnF4 phase), that the optimized Redlich- 
Kister coefficients for one and the same system can be 
different Q. 

Figure [3] shows the LiF-YF3 optimized phase dia- 
gram that was obtained in this work, together with the 
T - X experimental points for comparison. It may be 
noticed that the assessed values of the eutectic reaction 
temperature and the melting temperature of LiYF4 are in 
good agreement with experimental data. The assessed 
enthalpy and entropy at 298.15 K temperature are listed 
in table Q] The excess parameters assessed, given by Lq 
and L\ in the Redlich-Kister model (equation [5]l were 
L = -15883.80-24.187/ andLi =42271.89-23.267. 
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Figure 3: Theoretical phase diagram of the system LiF-YF3 calcu- 
lated assuming the experimental data from the literature (squares) and 
those measured in this work (stars). 



5. Conclusions 

Thermodynamic assessment has been performed on 
the L1F-YF3 binary system. The excess Gibbs pa- 
rameters for the liquid phase could be properly opti- 
mized using a Redlich-Kister polynomial model. The 
resulting theoretical phase diagram shows satisfactory 
agreement compared to the experimental phase dia- 
gram. Aff(298.15K) and S (298. 15 K) have been as- 
sessed for LiYF4, and a re-evaluation for AHf and Cp 
data for this intermediate compound was performed. It 
should be noted that the liquidus from the YF3 rich 



side meets the LiYF4 liquidus at a molar fraction of 
50% YF3. Results presented in this paper are contribu- 
tions to the more complete thermodynamic description 
of L1F-YF3 system and may be useful for the further 
thermodynamic assessment, in particular for ternary and 
multi component systems based on this binary phase di- 
agram. 
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